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We propose an experiment on quantum feedback control 
of a solid-state qubit, which is almost within the reach of the 
present-day technology. Similar to the earlier proposal, the 
feedback loop is used to maintain the coherent (Rabi) oscil- 
lations in a qubit for an arbitrary long time; however, this 
is done in a significantly simpler way, which requires much 
smaller bandwidth of the control circuitry. The main idea is 
to use the quadrature components of the noisy detector cur- 
rent to monitor approximately the phase of qubit oscillations. 
The price for simplicity is a less-than-ideal operation: the fi- 
delity is limited by about 95%. The feedback loop operation 
can be experimentally verified by appearance of a positive 
in-phase component of the detector current relative to an ex- 
ternal oscillating signal used for synchronization. 



The needs of quantum computing 1 are fueling a rapid 
progress in experiments with solid-state qubits. In par- 
ticular, quantum coherent (Rabi) oscillations have been 
demonstrated using superconducting charge, flux, and 
phase qubits 2 as well as double-quantum-dot qubits. 3 
Successful experiments with two superconducting qubits 
have also been demonstrated. 4 Even though at present 
only very basic operations with qubits are experimen- 
tally accessible, more advanced experiments are a natu- 
ral next stage. One of the directions for the advanced 
qubit control is realization of the quantum feedback con- 
trol of a solid-state qubit, 5 which can be used in a quan- 
tum computer for qubit initialization and is also an 
important demonstration by itself, clarifying the con- 
troversial issue of gradual collapse of a quantum state. 
(In optics quantum feedback control was proposed more 
than a decade ago 6 and has been already demonstrated 
exp er iment ally. 7 ) 

For the analysis of a quantum feedback we have to 
take into account the process of continuous qubit col- 
lapse. Therefore, the conventional approach to con- 
tinuous quantum measurement 8 ' 9 is inapplicable, and 
it is necessary to use the recently developed Bayesian 
approach 10 or the equivalent (though much different 
technically) approach of quantum trajectories. 11 The 
possibility of a quantum feedback is based on the fact 
that measurement by an ideal solid-state detector (with 
100% quantum efficiency rf) does not decohere a single 
qubit, 10 even though it decoheres an ensemble of qubits 
because each qubit evolves in a different way. The ran- 
dom evolution of a qubit in the process of measurement 
can be monitored using the noisy detector output, with 
accuracy depending on 77, so that for an ideal detector 



(j] = 1) even the monitoring of qubit wavefunction is 
possible. An example of theoretically ideal solid state de- 
tector is 12 ' 10 the quantum point contact (77 comparable to 
1 has been demonstrated experimentally 13 ). The single- 
electron transistor is significantly nonideal 9 ' 10 ' 14 (rj <C 1) 
in the semiclassical "orthodox" mode of operation; 15 
however, it can reach ideality in some modes based on 
cotunneling or Cooper pair tunneling. 16 

Monitoring of the quantum state in real time can nat- 
urally be used for continuous feedback control of a quan- 
tum system. In the experiment proposed in Ref. 5 the 
quantum feedback is used to maintain quantum coher- 
ent (Rabi 17 ) oscillations in a qubit for an arbitrary long 
time, synchronizing them with an external classical sig- 
nal. This is done by measuring the noisy current I(t) in a 
weakly coupled detector and using the quantum Bayesian 
equations 10 to translate information contained in I(t) 
into the evolution of qubit density matrix p(t). After that 
p(t) is compared with the desired quantum state Pd{t), 
and the calculated difference is used to control the qubit 
Hamiltonian in order to decrease the difference. Notice 
that the measurement backaction necessarily shifts the 
phase of Rabi oscillations in a random way (adding it to 
dephasing due to environment); however, the information 
contained in I(t) is sufficient to monitor this change and 
therefore restore the desired phase. 

An important difficulty in such experiment is a neces- 
sity to solve the Bayesian equations in real time. More- 
over, the bandwidth of the line delivering I(t) to the 
circuit solving the Bayesian equations, should be signifi- 
cantly wider than the Rabi frequency f2 (otherwise the in- 
formation contained in the noise is lost). Unfortunately, 
these conditions are unrealistic for the present-day exper- 
iments with solid-state qubits. (The "direct" feedback 
also analyzed in Ref. 5 does not require solving Bayesian 
equations, but requires a wide bandwidth of the whole 
feedback loop.) 

In this paper we propose and analyze a much simpler 
way (Fig. 1) of processing the information carried by the 
detector current I(t). The idea is to use the fact that be- 
sides noise, I(t) contains an oscillating contribution due 
to Rabi oscillations in the measured qubit. Therefore, if 
we apply I(t) to a simple tank circuit (which is in reso- 
nance with f2), then the phase of the tank circuit oscil- 
lations will depend on the phase of Rabi oscillations. In- 
stead of using the tank circuit, almost equivalent theoret- 
ically procedure is to mix I(t) with the signal from a local 
oscillator (Fig. 1) in order to determine two quadrature 
amplitudes of I(t) at frequency f2, which will carry infor- 
mation on the phase of Rabi oscillations. Since diffusion 
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FIG. 1. Schematic of the proposed quantum feedback loop. 
Two quadrature components of the detector current I(t) are 
used to monitor approximately the phase difference between 
Rabi oscillations and a local oscillator, which is used to control 
the qubit parameter H . The phase can also be monitored 
using a tank circuit. Positive average in-phase quadrature (X) 
is an experimental indication of quantum feedback operation. 



of the Rabi phase is a slow process (assuming weak cou- 
pling to the detector and environment), the further cir- 
cuitry can be relatively slow, limited by the qubit dephas- 
ing rate, but not limited by much higher Rabi frequency. 
The simplicity of the information processing and rela- 
tively small required bandwidth are the main advantages 
of this proposal in comparison with Rcf. 5 The experi- 
ment can be realized using either superconducting 2 ' 4 ' 18 
or GaAs 3,13 technology. 

The idea of this proposal partially stems from the fact 
that in absence of feedback the qubit Rabi oscillations 
lead to a noticeable peak in the spectral density Si(w) of 
the detector current at ui s=s fi, with the peak-to-pedcstal 
ratio up to 4 times 19 (somewhat similar experiments have 
been reported recently 20 ). Since 4 is not a big number, 
one would expect quite inaccurate phase information car- 
ried by current quadratures and therefore poor operation 
of the feedback. Surprisingly, the quantum feedback op- 
erates much better than it would be expected from clas- 
sical analysis. 

Let us consider a "charge" qubit (either double quan- 
tum dot or single Cooper pair box) with Hamiltonian 

7~{-qb = (£/2)(4 c 2 - c\ci) + H(c[c 2 + 4 c i)> where c\ 2 
and Ci t 2 are the creation and annihilation operators in 
the basis of "localized" (charge) states, e is their energy 
asymmetry, and the tunneling H = H + Hfb(t) can be 
controlled by the feedback loop (-ff/b)- We assume the 
standard coupling 21,10 ' 19 between the charge qubit and 
the detector (quantum point contact or single-electron 
transistor). Instead of writing Hamiltonians explicitly, 
we will characterize the measurement by two levels of 
the average detector current, ii and I 2 , corresponding to 
the two charge states, by the detector output noise Si, 
and by the total ensemble-averaged qubit dephasing rate 
r due to detector backaction and environment. Assum- 
ing sufficiently large detector voltage and quasicontinu- 
ous detector current J(t), we describe the qubit evolution 
by the Bayesian equations 10 (in Stratonovich form) 

pu - -2Hlmp 12 + 2p u p 22 [I(t) - I Q ] AJ/Sj, (I) 
Pi2 = ispi2 + iH(p n - p 22 ) - 7P12 

-(pii-P22)pi2[/(t)-/o]AI/Sj, (2) 



where Ti — 1, AI = I\ — I 2 , Iq = (I\ + I 2 )/2, and 7 
I — (A/) 2 /4S7. The qubit dccohcrcncc rate 7 = jd + 7e 
is due to detector nonideality, jd = (ry -1 — 1)(AI) 2 /Si, 
and due to additional coupling with environment (7,=). 
The current J(i) = I + (pu - p 22 )AI/2 + has the 
noise component £(t) with the flat (white) spectral den- 
sity Si. Notice that in the case e = (which is assumed 
unless mentioned otherwise), we can disregard the evolu- 
tion of Re/?i2 (it becomes zero at t 3> T -1 ), so only two 
degrees of freedom are left, which may be parameterized 
as pu — p 22 = Pcos(fii + 0) and 2lmpi 2 = Psin(fii + <f>), 
where the feedback-maintained frequency fi (see below) 
is assumed to be equal (unless stated otherwise) to the 
bare Rabi frequency fi = (4-ffg + s 2 ) 1 ^ 2 = 2H . More- 
over, in the ideal case 7 = the state eventually becomes 
pure, 10 so that P = 1 and the evolution can be described 
by only one parameter <p(t). 

We assume that two quadrature components of the de- 
tector current (Fig. 1) are determined as 

X(t) = f [I(t') - I ] co&(M)e- ( - t - t '^ T dt l , (3) 

J — CO 

Y(t) = f [I(t') - I ] sin(nt')e-(*-*')/ T dt' , (4) 

J — 00 

where f2 is the local oscillator frequency applied to the 
mixer, and t is the averaging time constant. Similar 
formulas are also applicable to the case of a tank cir- 
cuit with the resonant frequency f2 and quality factor 
Q = fir/2. If the detector current would be a har- 
monic signal I(t) = I + P(AI/2) cos(fit + <po), then 
<t>o = — arctan((Y)/(AT)), so it is natural to use 



= - arctan(F/A:) 



(5) 



as a monitored estimate of the phase shift <j>(t) between 
the Rabi oscillations and the local oscillator ((. . .) means 
averaging over time) . 

Let us assume 7 = and analyze first how close is the 
estimate <j> m (t) to the actual phase <j>(t) without feed- 
back, in which case <f> evolves in a diffusive manner due 
to detector backaction. Figure 2 shows the rms phase dif- 
ference Acj) rms — {{<t> m — 4>) 2 ) 1 ^ 2 (solid lines) as a function 
of t for several values of the dimensionless qubit-detector 
coupling C = (A/) 2 / SiHq, calculated numerically using 
Monte Carlo simulation of the measurement process. 10 
At weak coupling, C < 1, the curves practically coin- 
cide, and the minimum A<p rms w 0.44 is achieved at 
r w 4S7/(A/) 2 = 1/r, which is expectable since T deter- 
mines the phase diffusion: 10 ' 11 ' 19 ([^(i) - (f)(0)} 2 ) /t = Y. 
At larger r, <f> m includes too much of irrelevant informa- 
tion from distant past, while at smaller r the quadrature 
amplitudes suffer too much from the noise. At r — > (as 
well as at r — > 00) A(f> rms — > tt/VS ~ 1.81 that corre- 
sponds to the uniform distribution of A(j> — cj> m — <j> (com- 
plete absence of correlation between 4> and <j> m ) within ±7r 
interval (all phases are defined modulo 27r). 
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FIG. 2. Dependence of monitoring inaccuracy A<f) rms on 
averaging time r without feedback for several values of cou- 
pling C. Dashed and dotted lines correspond to classical sig- 
nals. Inset: distribution of Afi for several r at weak coupling. 

It it important to notice that the calculated A<fi rms 
is significantly smaller than for a naive classical case, in 
which the noise is not correlated with diffusive evo- 
lution of <fi. The dotted line in Fig. 2 shows the result 
for such a case at weak coupling, which has a minimum 
A(j) rrns w 1.0 [actually, for this curve we even increased 
the signal, assuming I{t) - J = v^(AJ/2) cos{Q.t + cfi) + 
which corresponds to correct spectrum 19 ]. Even 
more surprisingly, at r > 2.5Si/{AI) 2 the inaccuracy 
A(f> rms m the quantum case is smaller than for the clas- 
sical noiseless case, £(t) = (dashed line), which means 
that the noise improves the monitoring accuracy. This 
quantum behavior can be understood from the phase evo- 
lution equation 5,10 

4> = ~[i{t) - i Q ] sm(nt + <i>)(ai/Si) + n -n. (6) 

This equation 22 shows that the quadrature component 
of the noise which shifts the observed phase (j> m , also 
shifts the actual phase cfi in the same direction. In other 
words, when the noise imitates oscillations, it forces the 
real Rabi oscillations to evolve closer to what is observed. 

Inset in Fig. 2 shows the distribution of A<fi in the 
weak coupling limit for several values of r. The distri- 
butions are significantly non-Gaussian with the central 
part significantly narrower then A<fi rms . It is interest- 
ing that the value r = ASi/(AI) 2 corresponding to the 
minimum A<fi rms , does not provide the highest peak of 
the Aifi distribution. To find the best r in this respect, 
let us compare the monitored phase evolution <fi m — 
-[I(t)-I }sm(nt + <fi m )/(X 2 + Y 2 ) 1 / 2 , withEq. (6). We 
would expect the best approximation of <fi by <fi m when 
(X 2 +Y 2 ) = (5//A/) 2 . Using the definitions (3)-(4) and 
the current-current correlation function 19 (I(0)I(t)) = 
(Si/2)S(t) + (A//2) 2 cos(m)exp[-(A/) 2 t/85 , 7 ], we ob- 
tain (X 2 + Y 2 ) = S7t[1/4 + 1/(1 + 8S//(A/) 2 t)] at fir > 
1 and C « 1, so the condition (X 2 + Y 2 ) = (5//A/) 2 
is satisfied at r(A/) 2 /5/ = (2/5)(\/41 - 1) » 2.16. This 
indeed corresponds to the largest peak of A<fi distribution 
(see inset in Fig. 2). 
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FIG. 3. Dependence of the synchronization degree D on 
the feedback factor F in ideal case (7 = 0) for several r. Ex- 
perimentally D can be measured via average in-phase current 
quadrature (X). Dashed line is for a classical feedback. 

Reasonably small difference between <fi and <fi m in ab- 
sence of feedback implies that we can expect decent op- 
eration of the quantum feedback loop in which the phase 
estimate <fi m is used for determining the feedback action. 
Similar to Ref. 5 we consider the feedback loop, which aim 
is to suppress the fluctuations of the Rabi phase, so that 
the goal is <fi(t) = (or as small as possible). It has been 
shown that this goal can be fully reached using the linear 
feedback rule Hfb(t)/Ho = —Fcfi{t), which requires exact 
monitoring of (fi; here we analyze the operation of the 
feedback loop with Hf^/Ho = —F(fi m {t), where F is the 
dimensionless feedback factor (by definition \<fi m \ < 7r). 

We will characterize the performance of the feedback 
loop by the synchronization degree D — {Pit) cos <fi{t)) = 
IT — 1 where T = {Tr p{t)pd{t)) is fidelity and pa corre- 
sponds to the desired perfect Rabi oscillations {Pd = 1, 
(fid = 0)- Figure 3 shows (solid lines) the dependence of 
D on the feedback factor F for several time constants 
t in the case of weak coupling C = 0.1 and 7 = (we 
normalize F hy C, so the results practically do not de- 
pend on C for C < l). 23 One can see that each curve has 
a maximum, so that the "oversteering" effect at larger 
F makes the feedback performance worse (this is in con- 
trast to the case of Ref. 5 in which larger F was always 
better). Somewhat unexpectedly, r = 1/T = 4Si/{AI) 2 
is no longer an optimum, and smaller time constants are 
actually better. It can be shown that the feedback loop 
can operate even at r -C £1 _1 <C T _1 ; however, we are 
not interested in this regime because it requires a wide 
bandwidth of the control circuitry. Limiting ourselves 
to t ~ Si/{AI) 2 , we see that the maximum achievable 
synchronization degree D max is about 90% (that corre- 
sponds to the fidelity T of about 95%). It is impossible 
to reach 100% because the monitored simple phase esti- 
mate (fi m is significantly different from the actual (fi. It 
is interesting to note that a very crude estimate of D max 
as cos{A(fi rrns ) using min(A0 rms ) ~ 0.44 from the anal- 
ysis without feedback, works quite well, cos(0.44) = 0.90 
(though for different r). Dashed line in Fig. 3 shows the 
feedback performance for a classical signal corresponding 
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FIG. 4. Solid lines: synchronization degree D (and 
in- phase current quadrature (X)) as functions of F for sev- 
eral values of the detection efficiency n e f / . Dashed and dotted 
lines illustrate the effects of the energy mismatch (e ^ 0) and 
the frequency mismatch (fl / Oo)- 

to the dotted line in Fig. 2, assuming t(AI) 2 /Si = 1. As 
expected, it operates much worse than the quantum feed- 
back because of the reason discussed above. [The crude 
estimate D max ~ cos(A0 rms m i n ) = cos(l.O) = 0.54 still 
works well.] 

An important question is how the operation of the 
quantum feedback loop can be verified experimentally. 
One of the easiest ways is to check that the average value 
(X) of the in-phase quadrature component X(t) becomes 
positive, while in absence of feedback (F = 0) positive 
and negative values of X are obviously equally probable. 
Notice that any Hamiltonian control of a qubit which is 
not based on the information obtained from the detector 
(i.e. feedback control) cannot provide nonzero (A). 24 It 
is easy to show that (A) = [L> + (Pcos(2fi m + 0))]rA//4, 
and since the second term in brackets vanishes at weak 
coupling (and e = 0), therefore (A) is directly related 
to D. The numerical results for (A)/(tAJ/4) practically 
coincide with the curves for D in Fig. 3 (within the thick- 
ness of the line). 

The ideal case 7 = is obviously not realizable in 
an experiment because of finite nonideality of a de- 
tector (77 < 1) and presence of an extra environment 
(j e > 0). Both effects can be taken into account si- 
multaneously introducing effective efficiency of quantum 
detection r] e ff = [r/^ 1 + -f e Si/(AI) 2 ]~ 1 . Figure 4 shows 
(solid lines) the feedback performance for several values 
of i] e fj assuming t(AI) 2 /Si = 1. One can see that 
Veff ~ 0.1 is still a sufficient value for a noticeable oper- 
ation of the quantum feedback loop. Notice that D max 
is limited at least by the state purity, D max < P, which 
is 22 P w y/2rj e f f at r] eff < 1 and C/rj < 1 (D max = P 
can be reached by the feedback of Rcf. 5 but not by the 
feedback studied here). 

Finally, let us discuss how accurately the conditions 
O = Oo and e = should be satisfied in an experiment. If 
O is different from O , then without feedback the phase 4> 
linearly grows in time [Eq. (6)]. However, if the feedback 
loop operation is faster than |AO| = |0 — Oq|, the linear 



growth of <f> is stopped by adjusting the Rabi frequency 
O to match the desired frequency O. Dotted lines in Fig. 
4 show the feedback operation for r/ e ff = 0.2 and two 
values of AO, confirming still good operation at | AO| -C 
CO <~ T <~ r _1 . Notice that the frequency mismatch 
leads to nonzero {(j> m ) and therefore can be noticed and 
corrected. Energy mismatch (e 7^ 0) also worsens the 
performance of the feedback loop; however, the dashed 
lines in Fig. 4 (rj e ff = 0.5) show that a relatively large 
mismatch (e < H ) can be tolerated. 

In conclusion, we have proposed and analyzed the 
quantum feedback loop for a solid-state qubit, based on 
monitoring the phase of Rabi oscillations via quadrature 
components of the current in a weakly coupled detec- 
tor. Surprisingly, it operates much better than one could 
guess from relatively small spectral peak of the detector 
output in absence of feedback. The work was supported 
by NSA and ARDA under ARO grant DAAD19-01-1- 
0491. 
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